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Abstract 

Using an interface displacement model we calculate the shapes of nanometer- 
size liquid droplets on homogeneous cylindrical surfaces. We determine effec- 
tive contact angles and line tensions, the latter defined as excess free energies 
per unit length associated with the two contact lines at the ends of the droplet. 
The dependences of these quantities on the cylinder radius and on the volume 
of the droplets are analyzed. 
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I. INTRODUCTION 



The wetting properties of a fiber in liquid matrices (e.g., dye mixtures, polymer melts, or 
molten resins) play an important role in the textile industry and in the fabrication of high- 
performance, fiber-reinforced composite materials. Since contact angles of liquid droplets 
on solid substrates provide a valuable characterization of such wetting properties there are 
numerous experimental and theoretical studies of the shape and the spreading of droplets 
deposited on a cylindrical substrate (see, e.g., Refs. PHT5|). The morphology of liquid 
drops on a fiber is particularly interesting insofar as on a planar substrate there is only 
one, spherical caplike droplet shape, whereas on a cylindrical substrate droplets may exhibit 
two, topologically different shapes, a "clamshell"- and a "barrel" -type one, depending on 
the droplet volume, the contact angle, and the cylinder radius PH3L The aforementioned 
studies deal with thick fibers and large drops, i.e., the length scales are /zm and larger. In 
this range the fluid structures are determined by macroscopic properties alone, i.e., volume 
of the liquid, surface tension a of the liquid vapor interface, Young's contact angle 9^, and 
radius R of the cylinder. 

However, with the discovery of nanotubes the interest in such fluid structures has shifted 
to much smaller scales. There are several applications for which these small solid-fluid struc- 
tures are very important, (i) For fabricating valuable composite materials involving nano- 
tubes their wetting by the liquid host matrix is necessary to couple the inherent strength of 
the nanotubes to the matrix, reinforcing materials or fillers for plastics and ceramics (JHJ . (ii) 
Nanotubes can be used as supports for heterogeneous catalysis or as templates for creating 
small wires or tubular structures by coating them with metals or metal oxides in the liquid 
state [17] or by attaching inorganic and organic moieties to the nanotube surfaces [18]. (iii) 



In order to use nanotubes as "nanostraws" potential candidates for exploiting such capil- 
larity must be screened by first seeing if the liquid wets the outside of nanotubes ]0| . The 
performance of the nanotubes as catalysts, adsorbants, and deodorants can vary depend- 
ing on whether they are composed of carbon, boron nitride, or oxides (SiC>2, AI2O3, V2O5, 



M0O3, T1O2) ||20|| . This variety demonstrates, that the substrate potential of these tubes can 
be regarded as a tunable parameter, (iv) By using nanotubes as nanotweezers it might 
be possible to grab and manipulate small liquid drops. For this application the substrate 
must be nonwettable. 

These small scales are comparable with the range of the substrate potential of the cylin- 
ders and of the molecular forces between the fluid particles adsorbing on them. Thus the 
droplets form under the action of the so-called effective interface potential u which accounts 
for the net effect of the competition between the forces among the fluid particles and the 



substrate potential [G2] . Accordingly the calculation of the corresponding deformed droplet 
shapes requires a more detailed theoretical description which takes the effective interface 
potential into account. To our knowledge there is only one, recent publication in which 
this effect of u on the droplet shape on fibers has been analyzed f2~3fl . It is the purpose 
of our study here to refine and to extend this analysis in various directions. If the radius 
R of the fiber reduces to a few nm, as it is the case for nanotubes, the effective interface 
potential itself will depend on R and thus deviate from that of the corresponding semi- 
infinite planar substrate used in Ref. [p^] . Accordingly we present a systematic analysis of 
the dependence of the shape of the droplets and their suitably defined contact angles on 
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both R and the droplet volume. This enables us to describe systematically the crossover in 
shape and contact angle between those of droplets on a cylinder and on the limiting case 
R — > oo of a planar substrate. We remark on how the structure of the effective interface 
potential, depending on whether it leads to first-order or continuous wetting transitions, 
influences the morphology of the droplets. We confine our analysis to barrel-type droplets 
and estimate their metastability against roll-up to the clamshell configuration. Finally we 
study two types of line tensions. The first one concerns the line tension of three-phase 
contact between liquid, vapor, and substrate emerging at the ends of macroscopicly large 
drops on fibers which reduces to the familiar line tension of the straight three-phase contact 
line on a planar substrate. The second excess free energy concerns the effective line tension 
associated with the circular shape of the three-phase contact line on a planar substrate as 
function of the droplet volume. These results are relevant for understanding how to extract 
line tensions from contact angle measurements. 

We are encouraged to present our refined analyses by recent experimental advances 
to determine droplet shapes such as microscopic interferometry [p4| |, ellipsometric mi- 



croscopy ||25|| , scanning polarization force microscopy [E6|, and tapping- mode scanning force 



microscopy [27]. These techniques allow one to resolve drop profiles on the submicrome- 



ter scale [^,0 down to the nanometer scale [57], both vertically and laterally. In view 
of the numerous important applications mentioned above it would be rather rewarding to 
extend the application of these techniques to nonplanar substrate geometries in order to re- 
solve experimentally the shape of droplets on fibers and tubes as presented in the following 
sections. 



II. THEORY 

A. Free energy functional 

In cylindrical coordinates the droplet surface is described by a function h(z) or l(z) of 
the coordinate z along the symmetry axis of the cylinder (Fig. 0(a)). We define h(z) and l(z) 
such that h(z) is the local separation between the liquid-vapor interface and the symmetry 
axis of the cylinder and l(z) = h(z) — R is the local separation between the cylinder surface 
and the liquid- vapor interface, i.e., the liquid layer thickness. The droplet is also symmetric 
with respect to a reflection at the plane z = 0. For large values of |z|, i.e., at large distances 
from the droplet center at z = 0, the liquid forms a thin wetting layer of thickness Iq = h — R 
around the cylinder. For reasons of simplicity h{z) is henceforth assumed to be a unique 
function of z, i.e., we do not consider contact angles 9 > 90°. The shape of the liquid-vapor 
interface enclosing the droplet is determined by the interplay of three physical quantitites: 
the Laplace pressure [Q 2aH generated by the mean curvature H of the interface with 
surface tension a, the capillary pressure induced by the finite droplet volume, and 



the disjoining pressure or, equivalently, the effective interface potential u c acting on the 
liquid- vapor interface ||22|| ; u c (l; R) is the cost in free energy per surface area to maintain 
a homogeneous wetting layer of prescribed thickness I covering the cylinder surface and 
can be expressed in terms of the underlying forces of the substrate and between the fluid 
particles In the absence of the effective interface potential, i.e., for large droplets the 
liquid- vapor interface is a minimal surface under the constraint of a prescribed volume, i.e., 
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it exhibits a constant mean curvature. The influence of the effective interface potential is 
most pronounced near the cylinder surface within the range of the substrate potential and 
leads to a deviation of the actual profile h(z) from the shape which is determined by the 
aforementioned constant mean curvature condition. On the other hand, in the limit of large 
separation from the cylinder surface the mean curvature is asymptotically constant because 
there the influence of the effective interface potential vanishes. 

Independent of the size of the droplet, for later purposes we define the "reference config- 
uration" (see Fig. |](a)) 

a ref (z) = a(z)&( Zl - \z\) + h &(\z\ - z x ) (2.1) 

where a(z) is that constant-mean-curvature surface which touches the surface h(z), h(z = 
0) = a(z = 0), and exhibits the same curvature H at the apex, i.e., the two principal radii 
of curvature 

(i + WQ))^ = i with h , * 

h"(0) h"(0) dz V ; 

and h(0) = R+l(0) of h(z) and, correspondingly, of a(z) (see Fig. |l](a)) at the apex position 
z = are identical. O denotes the Heaviside step function; ±z± are those values of z where 
a(z) and the homogeneous wetting layer h(z) = ho intersect. In this sense the values z = rkz\ 
define the positions of the two three-phase contact lines forming the ends of the droplet. 
The "apparent contact angle" 9 is defined by the intersection of the barrel-shaped part a(z) 
of the reference profile and the homogeneous wetting layer h(z) = ho (see Fig. |l](a)): 

9 = lim arctan(|a'(z)|). (2.3) 



This apparent contact angle 9 can be expressed in terms of the measurable quantities 
apex height 1(0) of the droplet, radius R of the cylinder, radius R + lo of the wetting film, 
and radius of curvature Ro (Eq. (j2.2|)) of the profile h(z) at the apex: 

R + l(0) (R + l(Q)f-(R + l )i ( I l\ 
R + lo 2(R + lo) \R + l(0) + Ro)' 1 ' 



Within an interface displacement model (see, e.g., Ref. |]30| ) the equilibrium interface 
configuration h(z) for a droplet of prescribed excess volume V ex minimizes the free energy 
functional 

F ex [h(z)]=F[h(z)]-F[h ] (2.5) 

oo 

= 2tt J dz (a (h(z) ^l + h'(z) 2 - ho) + R (u c (h(z) - R) - u c (h - R)jj 



under the constraint 



7T 



oo 

J dz (h 2 (z)-h 2 ) =V ex (2.6) 
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and the boundary conditions h(\z\ — > oo) = ho- We have defined F ex [h(z)] as an excess free 
energy with respect to the free energy F[h ] of the homogeneous wetting layer h(z) = h 
rendering a mathematically well-defined, finite expression. The first contribution to F ex is 
the excess free energy due to the increase of the liquid-vapor interface as compared with a 
homogeneous cylindrical shape. In general a itself depends on the curvature and thus on R 
(see, e.g., Sec. 2.2 in Ref. [^9| and references therein); in the following, however, we do not 



discuss explicitly this additional parametric dependence on R. The second contribution to 
F ex is the free energy generated by the effective interaction between the cylinder surface and 
the liquid-vapor interface, reduced by the corresponding free energy for the homogeneous 
wetting layer. Since the substrate is considered to be homogeneous, u c (l) depends only 
on the radial distance I = h — R from the substrate surface. The equilibrium separation 
ho = Iq + R of the homogeneous wetting layer from the cylinder axis minimizes the free 
energy F(h) = 27rL(Ruj c (h — R) + ah) where L 3> z\ is the macroscopic length of the 
cylinder. The constrained minimum of Eq. ( |2.5| ) is given by the unconstrained minimum of 
the surrogate functional 

T ex \h(z)\ = F ex [h{z)} + k^tc j dz (h 2 {z) - h 2 ) - V ex j . (2.7) 

The corresponding optimal profile h(z, k) renders the equilibrium profile h(z, V ex ) = 
h(z, n(V ex )) upon expressing the Lagrange multiplier k in terms of V ex by inserting h(z, k) 
into the left hand side of Eq. (|2.6| ) which yields the implicit relation V ex (K). In order to 
avoid a clumsy notation, in the following we denote h(z, V ex ) by h(z). The Euler-Lagrange 
equation corresponding to Eq. (|2.7|) reads 



h=h{z) 

(2.8) 



h"(z) \ n TT , , R du c (h-R) 

o tt^-, : ; - -, , , ( n ^ ln = 2aH(z) = -k - v ' 



h(z)(l + h'(z) 2 y/ 2 (1 + h'{z) 2 f/ 2 J w h(z) dh 



This equation describes the balance between the Laplace pressure on the left hand side 
and the capillary plus disjoining pressure on the right hand side. One has k < for any 
barrel-shaped droplet. 



B. Reference profiles 

The reference profile a(z) minimizes a similar surrogate functional: 

4*)] = (2« W v/TW+KVM) +const (2.9) 

-21 

with the constant independent of a(z) and the boundary conditions a(±zi) = ho- Equa- 
tion ( |2.9| ) follows from Eq. (|2.7| ) by omitting u c and replacing k by k*. The corresponding 
Euler-Lagrange equation describes the constant-mean-curvature surface given by 
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/ 1 a"(z) \ 

° -TT7 ttttttat - 7 0NQ/0 = 2aH(z) = -k*. (2.10) 

\a(z)(l + a'(z) 2 ) 1 /2 (i + a '(z) 2 ) 3 / 2 ) K ' K ' 

According to the definition of a(z) the Lagrange multiplier — k* has to be chosen such that 
this constant mean curvature of this surface equals the mean curvature at the apex of the 
actual surface h(z) = l(z) + R: 

- K - = -J_ + ± (2.11) 

a R + l(0) R V ; 

Rq (see Eq. ( |2.2|) ) and R + 1(0) are the principal radii of curvature at the apex of the 
actual surface h(z), determined by the former Lagrange multiplier n(V ex ). The solution of 
Eq. ( |2.10D is given implicitly by 



a ( z ) / o \ -1/2 



/*((?)'( 



C-y 



z + zx, \z\<z u (2.12) 



which fulfils the boundary condition a(—z\) = ho', this determines implicitly z\ in terms of ho, 
V ex , a, and u> c . The integration constant C is determined by a'(0) = due to the symmetry 
of a(z). The integral in Eq. (|2.12|) can be expressed in terms of elliptic integrals pT|p2| . 

When the drop is macroscopicly large (V ex = oo) it is appropriate to adopt a slightly 
different point of view. In this case not the center of the droplet but the position of one of the 
two three-phase contact lines, which are defined by the intersection of the asymptote a m (z) 
(the constant-mean-curvature surface appertaining to the macroscopic drop) and h(z) = ho, 
is fixed at z = (see Fig. |l](b)). The actual interface profile interpolates between, e.g., 
h(z — > — oo) = ho and h(z — > oo) = a m (z). This configuration describes a solid cylinder 
which is in contact with bulk vapor on the left hand side (z — > — oo) and with bulk liquid 
on the right hand side (z — > oo). The analysis of the internal structure of a three-phase 
contact line on a homogeneous, planar substrate is based on a similar configuration (see 



Refs. [gj] and |3J] and references therein). The interface profile diverges in the limit z — > oo: 
h(z — > oo) — > oo, but this divergence is not linear as in the case of the planar substrate. 
A macroscopicly large drop implies k — *■ 0. In this limit the volume constraint loses its 
meaning. Instead the state of the system is fixed by different lateral boundary conditions. 
In this case the solution of Eq. Q2.12 ) is given by 



a m (z) = a(z; k = 0) = C cosh f — — — J (2.13) 

with two integration constants C and D. a m (z) describes a rotational surface with minimal 
surface area. From Eq. ( |2.13| ) one can easily see that the divergence of the interface profile 
for macroscopic drops is exponential, a m (z —> oo) = ^exp((z — D)/C), rather than linear 
as on a planar substrate. The reference profile appertaining to the macroscopic drop is 

a r ef,m(z) = a m (z)Q(z) + h Q{-z). (2.14) 

The slopes at the intersection of the asymptote a m (z) and the homogeneous layer h(z) = h 
at z = defines the contact angle 9 m (R) = 9(R, V ex — > oo). R = oo corresponds to a planar 
substrate for which the interface profile diverges linearly in the limit z — > oo: 
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a m ,oo(z) - R = l + ztan9 OQ (2.15) 

with the macroscopic contact angle = 6 m (R — > 00) on the planar substrate. 9^ obeys 
Young's law cos = (a wg — a w i)/a where a wg and a w \ are the wall-gas and wall-liquid 
surface tensions, respectively; o wg — a w i = u c (l Q ; R = 00) is determined by the effective 
interface potential of the corresponding planar substrate (see, c.f., Subsec. [II Q ). On the 
cylindrical surface the contact angle 9 m (R) does not follow from similar thermodynamic 
considerations but follows from the numerical analysis of the full profile h(z) for large V ex 
(see, c.f., Sec. fTT|and Fig. g). 

The integration constants C and D in Eq. (|2.13|) can be determined from the conditions 
a m (z — 0) = R + lo and a'(z = 0) = tan# m so that 

a>m(z) — Rcos6 m cosh f — h arccosh —J . (2-16) 

\ R cos 9 m cos 9 m J 

The series expansion of this expression in terms of small z/Ris 

a m (z/R <^l)=R + l + z tan9 m (R) + 0(z 2 ). (2.17) 

In the limit R — > 00 the region where the higher order terms are relevant is shifted towards 
z = 00 such that, with 9 m (R — ► 00) = 9^, one recovers the linearly diverging asymptote 
a mi00 (z) (Eq. ( TO ). 



C. Effective interface potential 

For the same liquid layer thickness / the effective interface potential u c (l; R) of a cylin- 
der differs from that of a planar substrate u) p (l). The full expression uj c (1;R) is presented 
in Ref. [5l| as obtained from density functional theory and within a so-called sharp-kink 




approximation for the solid-liquid and the liquid-vapor interface profiles. For reasons of 
simplicity, here we use the leading order of a series expansion of u> c (l; R) in terms of d w /R 
where d w is the radial extension of the volume excluded for the fluid particles due to the 
repulsive part of the substrate potential: 

with h = I + R and 2-F1 hypergeometric functions. In the limit l/R — > one recovers the 
expression for the effective interface potential of the corresponding planar substrate: 

u c (l; R -> 00) =u c {l -> 0,R) =u p (l) = aT 2 + bl~ 3 + cT 8 . (2.19) 

However, the power-law decay of u c (l — > 00) for a fixed, finite cylinder radius R is 

u c (l ^oo;R) = ^a^ + G(r 4 ), (2.20) 
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i.e., one power faster than that for the corresponding planar substrate. 

At present there exists, to our knowledge, only one study concerned with the shapes of 
droplets on cylinders within the range of the effective interface potential between the cylinder 
surface and the liquid- vapor interface ]23|. However, in Ref. |23j the disjoining pressure 



n c (0 = —{R/ (R + 0) du c (l)/dl on the right hand side of the Euler-Lagrange equation 
as a whole rather than only the effective interface potential u c is replaced by the disjoining 
pressure of the corresponding planar substrate H p (l) = —du) p {l) / dl. In view of Eqs. ( |2 . 1 9| ) 
and (p.20| ), except for the factor R/(R + l), this corresponds to the short-distance expansion 
(l/R — ► 0) of the effective interface potential of the cylinder. This replacement of the 
disjoining pressure by that of the planar substrate is expected to yield numerically reliable 
results only for large cylinder radii and small liquid layer thicknesses / <C R. Therefore in 
Sec. |TTT| we test the quality of this approximation (as well as that of the replacement of uo c 
alone by cu p ). 

So far, due to the volume constraint, our considerations apply to nonvolatile liquids. 
For volatile liquids any droplet surrounded by a macroscopic reservoir of vapor phase is 
thermodynamically unstable against evaporation, leaving behind only the thin equilibrium 
wetting film. However, we expect that the actual nonequilibrium state of a condensating 
or evaporating liquid observed within a time scale that is small compared with the typical 
condensation or evaporation time can be described by solutions of Eq. ( |2.8j ) with V ex given 
by its momentary value. Only the interface configuration for k — 0, i.e., V ex = oo, which 
interpolates between a homogeneous wetting layer and an exponentially diverging profile, 
describes a bona fide thermodynamically stable state which can be maintained by imposing 
appropriate boundary conditions (see above) at liquid-vapor coexistence for the bulk fluid. 
The thermodynamic state, which in a grand canonical ensemble is defined by temperature 
and chemical potential, enters parametrically into the actual values of the effective interface 
potential uj c and the liquid-vapor surface tension a. 



III. SHAPES OF DROPLET SURFACES AND CONTACT ANGLES 



We solve the Euler-Lagrange equation fl2.8|) numerically for fixed values of k and for a 
given effective interface potential uj c (1); the value of k, in turn, determines the excess liquid 
volume V ex and allows us to establish the relation K(V ex ). As boundary conditions in the 
case k < (leading to droplets of finite size) we use that h(z) must approach the wetting 
layer thickness ho for large z and that h'(z = 0) = 0. The distance L/2, at which the system 
is cut off, is chosen large enough so that h(z = L/2) and h'(z = L/2) attain their asymptotic 
values ho and 0, respectively, within prescribed accuracy. The reference profile a re f(z) is 
then calculated numerically by solving the differential equation (|2.10| ) with k* determined 
by Eqs. Q2.21) and ( 2.10Q and with a(z = 0) = h(z = 0) and a'(z = 0) = 0, up to the point of 



intersection of a(z) and ho which defines the coordinate z\\ a re f(z > z\) = h$. The contact 
angle 9 is determined from Eq. ( [2.3|) and, as a crosscheck, from Eq. (|2.4j) . 

In all numerical calculations presented henceforth we set a = 3a s 2 , b = — 5<rs 3 , and 
c = 3cxs 8 such that s sets the length scale for the range of the effective interface potential 
(typically s w lnm). We divide both sides of Eq. ( |2.8| ) by a so that uj(l)/a is dimensionless 
and n/a has the dimension of an inverse length. Alternatively, instead of introducing s 
as above one can choose y/a/a as the basic length scale which describes the decay of the 
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effective interface potential; for our choice yaja ~ 1.73s. The effective interface potentials 
u; c (7; R) (Eq. fl2.18| )) and u) p (l) (Eq. fl2.19| )) for the above choice of coefficients are shown in 



Fig.g. 

As a first example we solve Eq. (|2.8|) with the effective interface potential uj p {1) of the cor- 
responding planar substrate (Eq. ( |2.19| )) and the potential coefficients given above. Figure |3| 
shows the profile of the droplet surface on a cylinder with radius R = 100s for us/a = —0.1. 
This choice of k leads to a small droplet with V ex ~ 1.46 x 10 4 s 3 (i.e., containing roughly 
10 7 fluid particles) whose liquid- vapor interface lies entirely within the range of the effective 
interface potential. Therefore the deviation of the profile from the asymptote a(z) extends 
up to the apex of the droplet. The model effective interface potential used here resembles 
a typical interface potential leading to first-order wetting on a planar substrate [P^j . The 
droplet surface crosses the reference profile and, upon approaching the apex of the droplet, 
it reaches the reference profile from below. 

Carroll @ has shown that, in the absence of the effective interface potential, the axisym- 
metric droplet configuration is only stable for 



R j cos6 - 3 [^§-) + 1>0 ' 

i.e., if the droplet is large compared with the diameter of the cylinder and if the contact 
angle is small. When the droplet volume decreases or the contact angle increases the ax- 
isymmetric droplet becomes metastable against a so-called "roll up" towards the "clamshell" 
configuration pE[| . Applying the stability criterion Eq. (|3.1|) to the interface profile shown 
in Fig. |3], we find that this barrel-type configuration is possibly metastable towards forming 
the "clamshell" shape. A definitive statement about the stability would require to refine the 
criterion in Eq. ( |3.1| ) by incorporating the effect of the effective interface potential. How- 
ever, the determination of the non-axisymmetric "clamshell" equilibrium shape requires a 
much larger numerical effort and is therefore beyond the scope of the present paper. One 
can define a critical value V eX)C such that for V ex > V eX)C the axisymmetric droplet is stable. 
Upon increasing R, V eXtC increases, too; V exfi — > oo in the limit R — ► oo. Only for V ex = oo, 
i.e., for macroscopic drops, and for contact angles smaller than 90° (as stated in Sec. [IT] here 
we do not consider the case 9 > 90°) the rotationally symmetric interface shape is stable for 
any value of R. 

Figure |] shows the droplet shape for the same choice of potential parameters and for 
the same cylinder with R = 100s, but for ns/a = —0.005. This choice of k leads to a much 
bigger droplet with V ex ~ 1.67 x 10 8 s 3 . The apex of the droplet is located at such a large 
distance from the cylinder surface that the effective interface potential is almost negligible. 
Therefore the application of Eq. (|3.1| ) is reliable; it yields that this particular droplet is 
indeed stable against "roll-up" . In the vicinity of the cylinder surface the absolute deviation 
of the interface profile from the asymptote is similar to that in Fig. ^. As compared with the 
situation shown in Fig. |3], the point where l(z) crosses the reference profile a(z) — R is shifted 
to the right and lies near the three-phase contact line at z — z\. For the model effective 
interface potential used here and in Fig. |3], in the region around the apex of the droplet the 
profile lies below the reference profile. These results are in accordance with the findings for 



the planar, homogeneous substrate |]33| , |34]| with the same type of interface potential. If, on 



the other hand, uj{1) corresponds to a system undergoing a continuous wetting transition, 
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i.e., exhibiting a single minimum without a potential barrier, the profile of the droplet shape 
approaches its asymptote from the outside without crossing it (see Fig. 8 in Ref. fl34|). 



Figure displays the effect of the replacement of the effective interface potential of a 
cylinder u> c (l] R) (Eq. ( |2.18|) ) by that of the corresponding planar substrate oj p {1) (Eq. (|2.19| )) 
for a droplet with us/a = —0.005 so that 1(0) « 361s and z\ ~ 453s on a thin cylinder 
with R = 20s. This droplet also satisfies the stability criterion in Eq. ( |3.1|) . For reasons 
of clarity in this figure we have plotted the difference Ah(z) between corresponding profiles 
instead of the profiles themselves. The influence of approximating u> c (l; R) by uj p (1) turns 
out to be rather small: the difference Ah(z) is at most of the order of s. It is even smaller 
in the case of smaller droplets and thicker cylinders: the quality of approximating u> c (l] R) 
by uo p (l) improves if the cylinder is thicker and the droplet is smaller. In Ref. ||2^] the 



disjoining pressure of a cylinder rather than the effective interface potential is replaced by 
its planar-substrate counterpart. This corresponds to replacing the term (R/h) du c /dl on 
the right hand side of Eq. ( |2.8| ) by du p /dl. The dashed line in Fig. [| shows the effect of 
this approximation on the surface profile. As expected, the quality of this approximation is 
worse than the substitution of the effective interface potential alone, although the difference 
Ah is still of the order of s. 

Figure |B| shows the apparent contact angles for the examples presented in Figs. |3] and f|, as 
well as for the same set of potential parameters but with R = 20s, R = 200s, and R = 500s 
as function of the excess liquid volume V ex . Upon increasing R the curves are shifted upwards 
and to the right. Due to the exponential divergence of the interface profile of a macroscopic 
drop (which is more pronounced for smaller R, see the discussion of Eqs. (|2.16 ) and (|2.17|) 
above and, c.f., Fig. §), the determination of the apparent contact angles for very large drops 
is, in particular for thin cylinders, numerically difficult. However, the data indicate that, 
for any value of R, 6{V ex ) is a monotonously decreasing function, with a vanishing slope 
d6/dV ex = at V ex = oo. The differences between the contact angles for different R are 
minimal at V ex = oo. For any R there is a sizeable increase of the apparent contact angle 
upon decreasing droplet size. 

Figure |7] shows the effect of replacing u p by uj c on the apparent contact angles for the 
system with R = 20s (compare Fig. ||). The difference between the contact angles calculated 
by using u c and u p is significant. It is much smaller for thicker cylinders whose contact angles 
are displayed in Fig. |[ However, the qualitative functional form of the dependence of 9{V ex ) 
is not affected by the replacement of u c by u p . 

For macroscopic drops (i.e., n = or, equivalently, V ex = oo) the Euler-Lagrange equation 
is solved with the initial value h(z = Li) = /i and a small initial slope h'(z = Li) (e.g., 
h'(z = Li) = 10~ 8 ); the initial slope h'(z = Li) = would yield the trivial solution h = h . 
In order to find the asymptote we determine the integration constants C and D in Eq. ( |2.13| ) 
such that a m (L 2 ) = h(L 2 ) and a^(L 2 ) = h'(L2) where z = L 2 is the coordinate up to which 
the differential equation is integrated numerically; the system size L 2 — L\ is chosen large 
enough so that upon further increase of the system size C and D remain unchanged within 
prescribed accuracy. The contact angle 9 m can be inferred from the value of C (Eq. 



Finally the coordinate system is shifted laterally such that the intersection of a m (z) and 
ho define the position z = (which corresponds to z = —Z\ for k ^ 0). As mentioned 
before, for macroscopic drops the rotationally symmetric configuration satisfies the stability 
criterion Eq. fl3.1p for any R. 
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The dependence of the liquid- vapor interface profiles of macroscopic drops on the cylinder 
radius R is shown in Fig. H using u p (l) and with the same set of interaction potential 
parameters as in the previous examples. In accordance with Eq. ( |2.16| ), the interface profiles 
for cylinders of finite thickness diverge exponentially in the limit z — > oo. In the limit 
R — > oo the region where higher-order corrections to the linear behavior (Eq. ( |2.17| )) are 
relevant is shifted towards z — ► oo such that in the limiting case R = oo corresponding to 
the planar substrate the linear divergence of the reference profile is recovered. Figure |](a) 
displays the apparent contact angles 9 m corresponding to the profiles shown in Fig. || as 
function of R. Upon increasing the cylinder radius R, 9 m approaches Young's contact angle 
#oo for the planar substrate as 9^ — 9 m ~ Br 1 . 

In Fig. g(b) we show the apparent contact angles of small droplets on a planar substrate, 
i.e., in the limit B — > oo but with V ex < oo. In this case the reference configuration is a 
spherical cap whose circular base has a radius z\ (see Fig. 0(a)). For large droplets cos# 
reaches Young's contact angle cos 9^ according to the Neumann-Boruvka equation 

y 

cos# — cos^oo = — (3.2) 

<JZ\ 

which allows one to determine experimentally the line tension of three-phase contact on 
a planar substrate by varying the droplet size. Figure |9](b) demonstrates that this linear 
relationship between cos# and zV 1 is valid only for z±/s > 500, i.e., for z\ > 500nm. From 
Fig. |9](b) one infers that cos 6* decreases more rapidly than predicted by Eq. ( |3.2|) . This 
behavior can be accounted for by an effective line tension T e ff(zx) which due to the circular 
bending of the three-phase contact line is larger than the value of the corresponding 
straight three-phase contact line. Similar results have been obtained by Dobbs ||37|| . 

9 and 9 m depend on the liquid-vapor surface tension a which in turn also exhibits a 
behavior a(B) — a(oo) ~ Br 1 . In accordance with the discussion in Subsec. [II A| , Fig. | does 
not yet take into account this indirect dependence of 9 and 9 m on B via a(B). 



IV. LINE TENSION 

As long as the size of a droplet is finite and fixed it is impossible to extract from the total 
free energy unambiguously and in a strict thermodynamic sense a line tension associated 
with the three-phase contact lines at the ends of the droplet because there are arbitrarily 
many ways to form the total free energy as a sum of various terms. However, well-defined 
line tensions emerge as coefficients in the size dependence of the free energy of droplets 
upon approaching macroscopic drops. To this end we consider the limit of large drops, i.e., 
V ex s~ 3 1 and Z\/s ^> 1 (see Fig. |TJ(a)). Within the interface displacement model the 
excess free energy in Eq. ( |2.5| ) can be rewritten as 

F ex = oA h - 4ttBuj c (Io)zi -aA c + C (4.1) 

where 

A b = 2n J dza ref (z)^l + (a' ref (z)) 2 (4.2) 

-21 
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is the surface area of the "barrel" part of the reference surface a re f(z), —Z\ < z < zx, 
A c = 4nh zi is the surface area of the cylinder with radius ho and length z x (see Fig. 0(a)), 
and 

C = 2ix I dz ( a 



h(z) y/l + {h'{z)f - a ref (z) yfl + (a' ref (z)y 



-zi 



For V ex — > oo, F ei . (Eq. ( |4.1| )) is dominated by the term aA b which scales proportional 
to the surface area of the drop and thus represents a two-dimensional contribution. The 
leading subdominant terms are — 4ttRuj c (1 )zi and —4irhoazi, which scale with the linear 
dimension 2zx of the drop representing one-dimensional contributions. Finally, the last term 
in Eq. ( |4.1| ) remains finite for V ex , A b , and Z\ — >• oo: 

C(z 1 )=T + 0(z^ 1 ) (4.4) 

and thus represents a zero-dimensional contribution. In its turn T depends on the cylinder 
radius R such that for large R it scales proportional to R which leads to the following 
definition of an excess free energy per unit length, henceforth called "line tension", associated 
with the two contact lines formed at the ends of the droplet with total length AtiR: 

7 (4.5) 



AixR 



We note that T is only well-defined in the thermodynamic limit lim [F ex — aA b + 

A(,,zi— »oo 

<AnRtjj c {lo)zi + crA c ] so that higher order terms, e.g., ~ zj~ , omitted on the right hand 
side of Eq. (|4.4|) drop out. In the following T is understood to have been obtained via this 
procedure. It can be expressed in terms of the solution of Eq. ( |2~8| ) for k = and of a re f >m (z) 
(Eqs. fln3D and (gig)): 



T = An J dz (a (h{z)y/l + (h'(z)f - a ref , m (z)Jl + {a' refm (z)f) 

— oo 

oo 

+ * - R) - -«)))+ *Rf«^M ~ «>. (4-6) 



Consequently, T is twice the characteristic excess free energy associated with the structure 
of a macroscopic drop near one of its ends without interference from the other end. On the 
other hand, T, and thus r, are defined for any value of R. The ratio r formed in Eq. ( |4.5| ) 
has the property that in the limit R — > oo it reduces to the line tension of the straight 
three-phase contact line on the corresponding planar substrate (see, e.g., Refs. [j^ , j^ , j34] ,j38|l ) , 
which is an experimentally observable quantity (compare Fig. |9](b)). 

At this stage one should note that the above considerations tacitly assume that another 
thermodynamic limit concerning the total system size, such as the volume of the surrounding 
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vapor phase and the length L of the solid cylinder, has already been carried out in advance: 
F(h Q ) is proportional to L and has been subtracted before. Moreover, we have not considered 
the bulk free energy of the surrounding vapor phase and the bulk free energy of the liquid 
in the drop proportional to V ex because they do not enter the description of the droplet 
shape in terms of an interface displacement model. As a careful analysis of the line tension 
^ within density functional theory for a volatile liquid at gas-liquid coexistence shows, 
in comparison with this more complete theory the interface displacement model misses a 
contribution which is independent of the shape l(x) and is determined by 9^ and 1$ (the 
first term in the sum in Eq. (2.19) in Ref. []34|] , denoted as f in Eqs. (4.2), (4.3), and (4.5) 
in Ref. ||33[| ). Since, however, this constant contribution turns out to be numerically much 
smaller than those contributions captured by the interface displacement model (see Fig. 15 
in Ref. P3"|), we have refrained from determining it for the present, much more complicated 



geometry, assuming that the size ratios of these types of contributions remain roughly the 
same for the planar and the cylindrical substrate. 

Figure |1(] shows the dependence of the line tension r on the cylinder radius R using 
the planar effective interface potential u> p (l). r is given by sa times a numerical factor of 
the order of 1. It decreases monotonously for decreasing R and attains its maximum value 
Too for R — ► oo as — t(R) ~ R~ l . We note that this decrease of the line tension upon 
decreasing the radius of curvature R of the contact line is opposite to the increase of the 
effective line tension T e ff{z\) observed for a decreasing radius of curvature Z\ of the circular 
three-phase contact line on a planar substrate as can be inferred from Fig. ||(b). Thus line 
tensions of curved three-phase contact lines can be smaller or larger than the line tension of 
the corresponding straight contact lines. 

Whereas r M is experimentally accessible by monitoring the apparent contact angle of 
sessile droplets on a planar substrate as function of the droplet size, r(i?) cannot be deter- 
mined experimentally by direct observation. The basic reason for this difference is that the 
length 2ttz\ of the three-phase contact line of the sessile drop on the planar substrate can 
vary as function of the droplet size so that the optimal shape of the droplet responds to 
the associated cost 2 r KZ\T 00 of the free energy, whereas the excess free energy 2itRt(R) for 
the ends of the droplets is a constant contribution with respect to the droplet size on the 
cylinder due to the fixed value of R. This, however, holds only for the "barrel" -type shape of 
the drop, for which the length of the three-phase contact lines is fixed. For "clamshell" -type 
droplet shapes the length of the three-phase contact line does depend on the droplet size 
so that in this case the line tension will influence the droplet shape. Nonetheless there are 
systems for which t(R) can be experimentally relevant. If the cylinder is not a hard solid 
rod but consists of a soft material like, e.g., vesicles or tobacco viruses which float vertically 
at the liquid- vapor interface of a solvent, the positive line tension t(R) will strangle the 
object locally, depending on its restoring elastic forces. According to Fig. |10| this tweaking 
force —d(T(R)R)/dR is weaker for thin cylinders. 



V. SUMMARY 



Based on an interface displacement model (Eq. ( ^.5[ )) we have analyzed the shape and 
the free energy of "barrel" -type droplets of fixed volume V ex covering a cylindrical substrate 
of radius R (Fig. [TJ). For sufficiently small droplets their shape h(z) = R + l(z) is not 
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only governed by the surface tension a of the liquid- vapor interface but also by the effective 
interface potential u c (l; R) with a generic form as shown in Fig. ||]. We have obtained the 
following main results: 

1. Figures |] and ^ show how the deviation of the actual droplet shape h(z) from a suitably 
defined reference configuration a re f(z) depends on the droplet size. The reference 
configuration is uniquely defined by the requirement to touch the actual shape at 
the apex and to have a constant mean curvature which equals the actual one at the 
apex. a re f(z) allows one to introduce an apparent contact angle 9, characterizing 
the actual shape, which can be expressed in terms of the experimentally accessible 
quantities cylinder radius R, radii of curvature at the apex, height 1(0) of the droplet, 
and thickness Iq of the wetting layer outside the barrel (Eq. ( |2.4|) ). 

2. The dependence of the effective interface potential uj c (1; R) on the cylinder radius R 
influences the shape of the droplet on the scale of the range s of u c (l; R) (Fig. [5|); this 
dependence has a rather marked effect on the apparent contact angle (Fig. |7|). 

3. The apparent contact angles increase for smaller droplets and for thicker cylinders 
(Fig. The contact angles of macroscopicly large drops approach Young's contact 
angle on a planar substrate proportional to 1/R (Fig. ||(a)). 

4. In the limiting case of small droplets on a planar substrate the circular bend of the 
three-phase contact line leads to an effectively increased value of the corresponding 
line tension (Fig. ||(b)). This deviation from the Neumann-Boruvka equation becomes 
relevant if the droplet radius is less than roughly 500nm. This observation is relevant 
for experimental determinations of line tensions via contact angle measurements. 

5. For macroscopicly large drops their shape l(z) increases linearly on a planar substrate 
but exponentially on a cylinder (Fig. 0(b)). Figure || illustrates the smooth crossover 
between these types of behavior for increasing cylinder radii. 

6. For large cylinder radii the line tension associated with the ends of macroscopicly 
large drops approaches the line tension of three-phase contact on a planar substrate 
proportional to 1/R (Fig. HIT). 
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FIGURES 




FIG. 1. (a) Schematic longitudinal cross-section through a liquid droplet residing on a homo- 
geneous cylindrical substrate with radius R. The system is rotationally symmetric around the z 
axis and symmetric with respect to a reflection at the plane z = 0. h(z) (full line) describes the 
liquid-vapor interface profile; the thickness of the liquid layer is l(z) = h(z) — R. The droplet shape 
is determined by the balance of the Laplace, the capillary, and the disjoining pressure (Eq. Q2.8D ). 
Far from the droplet center, i.e., in the limit \z\ — ► oo, the profile h(z) reduces to a homoge- 
neous layer h{\z\ — » oo) = ho- a re j(z) (dashed line) describes the reference surface which consists 
of a surface a re f(\z\ < z\) = a(z) with constant mean curvature and of the homogeneous layer 
GVe/(M > z i) = ho. a{z) is determined by a(z = 0) = h(z = 0) (i.e., it touches the actual surface 
profile h(z) at the apex) and by the condition that the two principal radii of curvature h(0) and 
Rq (which is the radius of curvature of the planar curve (z,h(z)), see Eq. ( |2.2[) ) of the actual 
and the reference surface at the apex are identical. The break in the slope of a re f(z) at z = ±z± 
defines the apparent contact angle 6. (b) Same as in (a), but for a macroscopic drop, i.e., infinite 
excess volume V ex . Choosing the position of one of the contact lines as the origin z = leads to 
a configuration for which the cylinder is in contact with bulk vapor for z — > — oo and with bulk 
liquid for z — > oo. a re f >m (z — > oo) and h(z — > oo) diverge exponentially. 9 m is defined by the break 
in the slope of a re f^ m {z) at z = 0. 
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FIG. 2. Model effective interface potentials u(l) in units of the planar liquid-vapor surface 
tension a as used in all numerical calculations. uj c (1;R) and u> p (l) are given by Eq. fl2l|) and 
( 2,19| ), respectively, with a/a = 3s 2 , b/a = —5s 3 , and c/a = 3s 8 where s sets the length scale of 
the system (typically s ~ lnm). The full line denotes the effective interface potential lo p {1) for the 
planar substrate, the dashed line denotes oo c {l]R) for a cylinder with R = 20s. For this choice of 
parameters and within the range of values of I shown here the effective interface potentials even of 
thin cylinders barely differ from that of a planar substrate. Only for large I the long-range decay of 
uj p {l) and to c (l) differ qualitatively (compare Eq. ( |2.20| )). This type of effective interface potential 
with a global minimum at Iq/ s ~ 1.3 and a second, local, minimum at / = oo leads to a first-order 
wetting transition of a planar substrate at a higher wetting transition temperature at which the 
first minimum raises up to w = 0. In the case of a continuous wetting transition would exhibit 
a single minimum and approach u = from below in the absence of a potential barrier in between. 
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FIG. 3. Profile l(z) of the droplet surface (full line) and the corresponding reference profile 
a re f(z) — R (dashed line) for R = 100s, ns/a = —0.1 so that V ex ~ 1.46 x 10 4 s 3 , and the effective 
interface potential u> p (l) of the corresponding planar substrate as shown in Fig. & The droplet is 
so small that its liquid-vapor interface lies entirely within the range of uj p (1). For z/s > 5.25 the 
interface profile l{z) lies above the reference profile a{z) — R and, for z/s < 5.25 upon approaching 
the apex of the droplet, it reaches the reference profile from below. This particular droplet is 
possibly metastable against the "roll up" to the "clamshell" shape because the stability criterion 
Eq. (3.1) is not satisfied. 
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FIG. 4. Profile l(z) of the droplet surface (full line) and the corresponding reference profile 
a re f(z) — R (dashed line) for the same set of parameters as in Fig. |||, but with ns/a = —0.005 so 
that V ex ~ 1.67 x 10 8 s 3 . In contrast to the situation shown in Fig. |3|, here the droplet is so large 
that it extends up to distances from the cylinder surface where the effective interface potential is 
negligible. The inset magnifies the region around the three-phase contact line at z = z\ ~ 483s. 
In the region around the apex of the droplet the profile l{z) lies below the reference profile. Since 
the effective interface potential is the same as in Fig. |3|the absolute deviation of h(z) from a re f(z) 
is about the same size as in Fig. |3] (see the inset). According to Eq. ( [O ) this droplet is stable 
against "roll up". 
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FIG. 5. Full line: difference Ah(z) = h c (z) — h p (z) between the two droplet profiles h c and 
hp for Ks/a = —0.005 (leading to a droplet with 1(0) ~ 361s and Z\ « 453s) and R = 20a which 
are calculated by using uj c (1; R = 20s) and u> p (l), respectively. For comparison, the dotted line 
shows the difference between the profile h c for R = 20s and ns/a = —0.005 and the profile h p for 
R = 20s but with k chosen such that the excess volumes of liquid for both profiles h p and h c are 
identical. In order to achieve this, for the determination of h p the value of ks/o has to be increased 
by 1.63 x 10~ 6 . In both cases the maximal difference is of the order of s. Since Ah changes sign 
as function of z, it is not possible to find an effective value K e ff such that for a given value of k 
the resulting profile calculated with ui p is identical with the profile h c . The dashed line denotes 
the difference Ah between the profile h c (z) for us/a = —0.005 and R = 20 and the profile h p (z) 
for which, as in Ref. [p3| , in Eq. the entire disjoining pressure (R/(R + I)) dw c /dl (instead 

of only duj c /dl) is replaced by the disjoining pressure of the planar substrate dw p /dl; here for the 
determination of h p the value of ns/a has been increased by 6.58 x 10 -6 in order to have identical 
excess liquid volumes. This approximation is worse than the substitution of the effective interface 
potential alone, although the difference Ah is still of the order of s. 



21 




FIG. 6. Apparent contact angles 9 for the same effective interface potential u p (l) as used for 
Figs. H] and ^| and for different cylinder radii R as function of the liquid excess volume V ex . The 
main effect of increasing R is a shift of the curves upwards and to the right. The limit V ex — > oo 
corresponds to macroscopic drops or, equivalently, a cylinder in contact with bulk vapor at one 
end and with bulk liquid at the other (compare Fig. |](b) and, c.f., Figs. || and ||). The differences 
between the curves are minimal in the limit V ex — > oo. 9(V ex ) is, for any value of R, a monotonously 
decreasing function with vanishing slope d9/dV ex = in the limit V ex — ► oo. The symbols ♦ indicate 
the critical values V ex , c above which according to the criterion in Eq. the axisymmetric droplet 
shape considered here is stable, but below which the "clamshell" configuration is stable. 
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FIG. 7. Apparent contact angles 9 calculated by using uj c (1; R) (dotted line) as compared with 
those calculated based on uj p (1) (full line) as function of the excess liquid volume V ex for a cylinder 
with R = 20s (i.e., the full line is identical with the full line in Fig. |6|). Both curves exhibit the 
same qualitative behavior. For thicker cylinders whose contact angles 0(V ex ) are shown in Fig. || 
the differences between the curves are much smaller. The symbols ♦ indicate the critical values 
V exc above which the axisymmetric droplet shape considered here is stable, but below which the 



"clamshell" configuration is stable (compare Eq. (3.1)). The difference between the two curves is 
largest for macroscopicly large drops. 
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FIG. 8. (a) Liquid-vapor interface profiles l{z) = h(z) — R of macroscopic drops (V ex = oo) for 
R = 20s (full line), R = 100s (dotted line), R = 1000s (dashed line) and R = oo (dashed-dotted 
line), calculated with to p (l). (b) The profiles l(z) (thick lines) and their corresponding reference 
configurations a re f(z) — R (thin lines) for R = 20s and R = oo are shown in the magnified region 
around the contact line at z = 0; here we use same graphical notation as in (a). Each reference 
configuration consists of the horizontal line l{z) = 1$ for z < and of the asymptotic branch 
dm(z) — R for z > 0; for R < oo this asymptote diverges exponentially in the limit z — > oo (see 
Eq. ( [2.16D ). R = oo corresponds to the planar substrate for which a mj00 (z) is a linear function. As 
demonstrated in (a), in the limit R — > oo the region where higher-order corrections to the linearly 
diverging asymptote become relevant is progressively shifted towards z — > oo, such that for R = oo 
only the linear divergence remains. The appertaining contact angles 9 m and the line tension r 
as function of R are shown in, c.f., Figs. [9| and [l(], respectively. Since for macroscopic drops the 
contact angles attain a finite value, for any radius R the stability criterion in Eq. ([O]) for the 
barrel-type shape is fulfilled. 
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FIG. 9. (a) Apparent contact angles 9 m for macroscopic drops (V ex = oo) as function of 
R (full line) corresponding to the interface profiles shown in Fig. ||. As indicated by the dashed 
line, in the limit R — > oo the apparent contact angles attain Young's contact angle #oo ~ 30.5° 
on the corresponding planar substrate as — m ~ R , (b) Apparent contact angle cos$ of 
small droplets (full line) on the planar substrate as function of the radius z\ of the base of the 
spherical cap acting as the corresponding reference configuration. For large drops cos 9 varies 
according to the Neumann-Boruvka equation (dashed line) cos 9 — cos#oo = —^/(azx) which 
allows one to determine experimentally the line tension of three-phase contact on a planar 
substrate. Here Too = 1-31<ts. The deviation between the full line and the dashed line shows that 
the Neumann-Boruvka equation is applicable only for z\/s > 500. 
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FIG. 10. Dependence of the line tension r (full curve) of macroscopic "barrel"-type drops 
(k = 0) on cylinders of radius R. The corresponding interface profiles and apparent contact angles 
are shown in Figs. || and 0, respectively, r is calculated by using the planar effective interface 
potential u> p (l). As indicated by the dashed line, r approaches Too for R — > oo as — r ~ R . 
Too is the line tension of the straight three-phase contact line on the corresponding planar substrate. 
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